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We analyze a sample of W-\- jet events collected with the Collider Detector at Fermilab (CDF) 
in pp collisions at ^/s = 1.8 TeV to study it production. We employ a simple kinematical variable 
Ti,, denned as the scalar sum of the transverse energies of the lepton, neutrino and jets. For events 
with a W boson and four or more jets, the shape of the Ti, distribution deviates by 3.8 standard 
deviations from that expected from known backgrounds to it production. However this distribution 
agrees well with a linear combination of background and it events, the agreement being best for a 
top mass of 180 GeV/c 2 . 

PACS numbers: 14. 65. Ha, 13.85.Ni, 13.85.Qk 



The existence of the top quark has recently been estab- 
lished by the CDF [1,2] and DO [3] collaborations. In the 
CDF analyses [1,2], b quark tagging was used to select 
it candidates in a sample of W+ > 3-jet events, where 
the W decays into ev or \iv. A W+ > 4-jet subsample 
was then used to reconstruct the top mass under the hy- 
pothesis that the top quark decays into a W boson and 
a b quark. The mass peak in the 6-tagged W+ jet data 
gives evidence for the top quark using a kinematic vari- 
able. In addition a study [4] of the transverse energy 
(Et) [5] spectrum of the second and third highest Et 
jets was able to identify a contribution of ti production 
in our data. 



For this analysis we use the CDF W+ > 4-jet sub- 
sample and study the variable Ti, which is defined as the 
scalar sum of the lepton Et, the neutrino Et as mea- 
sured by the missing- Et ($t) m the event, and the Et of 
each jet [6,7]. The variable Ti is strongly correlated with 
the center of mass energy (■s/i) of the parton-parton hard 
scattering process. In it events, it is also correlated to the 
transverse mass [8] of the it system. For a top mass larger 
than about 140 GeV/c 2 , Ti has discriminating power be- 
tween it and W + multijet background events, because 
the W and b quark resulting from the top quark decay 
have higher transverse momenta (Pt) than radiated glu- 
ons in background processes. For this study no 6-tagging 
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or top event reconstruction is required. Therefore, this 
analysis has larger acceptance for a it signal and has dif- 
ferent systematic uncertainties. It is not affected by am- 
biguities in jet-parton assignments, as is the case for the 
event-fitting algorithm used to determine the top mass in 
Ref. [2] and [3]. In addition, it is sensitive to unconven- 
tional top decay modes which do not involve b quarks, 
such as t — > Ws or t — » Wd. In these respects the study of 
T~L supplements Ref. [2] by providing additional informa- 
tion on the top quark mass and the it production rate. 

We report here on a study of the W+ > 4-jet data sam- 
ple used in Ref. [2] and based on 67 pb _1 of integrated lu- 
minosity. We search for it events where one W decays lep- 
tonically into a Iv pair (l = e, /x), and the other W decays 
into quarks. The CDF detector and the requirements on 
the lepton and^y are described in Refs. [1,2] and result 
in approximately 71,500 events. We then require events 
to have 3 jets that pass "high-threshold" cuts, with the 
observed calorimeter Et> 15 GeV and \rjj t t\ < 2.0, and 
at least one additional jet that passes a "low-threshold" 
cut of observed Et > 8 GeV and \r)j & t\ < 2.4. These jet 
cuts provide the 99-event signal sample used for this pa- 
per and for the mass analysis in Ref. [2], out of which 88 
events had a good reconstructed fit to the it hypothesis. 
When calculating 7i, the $t and the Et of the above 4 
jets are corrected for detector response. For 35% of the 
events in the W+ > 4-jet sample, there are extra jets 
that pass the low-threshold cuts and the corrected Et of 
these extra jets are added to the T~L variable. 

The dominant background to it production in the 
W+jets mode is direct production of a W which re- 
coils against energetic light quarks and gluons. Other 
backgrounds involving real W bosons are: WW, WZ +]ets 
where a W decays into lv\ and W+jet events where 
W-^tv and T^ti>iv T . The non-W backgrounds are: QCD 
multijet production where one jet fakes an electron or a 
muon (QCD fakes); bb + multijet production where one 
of the b quarks decays semileptonically; ZZ+jets with one 
Z decaying leptonically but with only one lepton found; 
Z — > tt followed by r— > lvtv T ; and Drell-Yan (j*, Z°) 
production of lepton pairs along with extra QCD jets. 

We model the shape of % for the W+ jets background 
using the VECBOS [9] Monte Carlo program. From the 
data and various Monte Carlo estimates, we obtain the Ji 
spectrum for all other background processes and find it 
is well matched to the VECBOS Ti. spectrum. The main 
systematic uncertainty in the VECBOS calculation is the 
Q 2 scale, which determines the strong coupling constant 
and therefore the production rate and the shape of the 
kinematic distributions. We use two choices of Q 2 scale, 
the square of the average Pt of the jets ((Pt) 2 ) and the 
square of the W boson mass (M^). The tt Monte Carlo 
samples are generated with the HERWIG program [10]. 
All the Monte Carlo samples are processed by a detector 
simulation program and reconstructed in the same way 
as the data. We emphasize that this analysis is based 
on comparing the shape of the measured T~L distribution, 
not the rate, with the signal and background predictions. 



TABLE I. Definition of the W+ jets control and signal 
samples. The third, fourth and fifth columns list the criteria 
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To test how well the VECBOS calculation models the 
background, we define two top-depleted W+3-jet control 
samples, and a low-threshold W+ > 4-jet signal sample 
in addition to the high-threshold signal sample defined 
above. These 4 samples are summarized in Table I. The 
it contamination in the control samples is expected to be 
about 1% (10%) for the low- (high-) threshold data [1], 
assuming a top mass of 175 GeV/c 2 . These two con- 
trol samples provide a means of verifying the W+ 3-jet 
background calculation. The low-threshold signal sample 
with T~L < 200 GeV is depleted in top events and allows 
us to test the W+ > 4-jet background calculation. 

We first compare the data and the VECBOS prediction 
for our control samples. The 7i plots for both control 
samples are shown in Fig. 1, along with the VECBOS 
calculations for Q 2 = (Pt) 2 mixed with the expected top 
component (1% top in Fig. la and 10% top in Fig. lb). 
The Monte Carlo distributions describe the data well. 
For Q 2 = M$y, we find that the low threshold VEC- 
BOS distribution is slightly shifted to larger values of 
T~L. However, the size of this shift is similar to the un- 
certainty in T~L from the 10% systematic error in the jet 
energy scale [1]. We therefore cannot reject the Q 2 = M^ 
scale on the basis of this comparison. The high-threshold 
data agree equally well with the VECBOS predictions 
for Q 2 = (P T ) 2 and Q 2 = M^. We conclude that the 
VECBOS calculation, plus the small it contamination, 




200 400 600 200 400 600 



H (GeV) H (GeV) 

FIG. 1. Comparison of Ti, distributions for the control 
data (solid line) and the VECBOS (Q 2 = (P T ) 2 ) prediction 
(shaded) for (a) W+ 3-jet events passing the low-threshold 
cuts, and (b) W+ 3-jet events passing the high-threshold cuts. 
The VECBOS prediction plus 1% top for (a) and 10% top for 
(b), has been normalized to the data. 



3 



CDF3110_v21 



FERMILAB-PUB-95/149-E 



June 13, 1995 




20 



15 



10- 



5- 



(b) 



— Data 
□ VECBOS 



46 



600 



200 400 600 200 400 

H (GeV) H (GeV) 

FIG. 2. Comparison of Ti, distributions for the signal data 
(solid line) and the VECBOS (Q 2 = (Pt> 2 ) Monte Carlo 
prediction (shaded) for (a) W+ > 4-jet events passing the 
low-threshold cuts, and (b) W+ > 4-jet events passing the 
high-threshold cuts. The Monte Carlo prediction is normal- 
ized to the expected number of background events obtained 
by the fitting procedure described in the text. 

describes both the low- and high-threshold control data 
well within our systematic uncertainties. 

We next compare the data and the VECBOS predic- 
tion for our signal samples. The Ti. plot for W+ > 4-jet 
events passing the low-threshold cuts is shown in Fig. 2a, 
together with the VECBOS prediction for Q 2 = (P T ) 2 . 
The VECBOS distribution is normalized to the expected 
number of events obtained by the fit described in the 
next paragraph. The data and VECBOS prediction agree 
well in both the peak position and the shape for Ti. be- 
low 200 GeV, suggesting that the low-7^ events are pre- 
dominantly background. However, on the high-7^ side 
above 200 GeV, a shoulder is seen in the data above 
the VECBOS curve. The histogram of Ti. for the 99 
W+ jet events passing the high-threshold cuts is shown 
in Fig. 2b. This sample consists of less than half the 
events in the low-threshold sample and is expected to 
include only 5% less top. The corresponding VECBOS 
distribution with Q 2 = (Pt) 2 is also shown, normalized 
in the same way as for the low-threshold sample. The 
data is significantly broader than the VECBOS predic- 
tion. Performing an unbinned Kolmogorov-Smirnov test, 
we find that the probability for the Q 2 = (P T ) 2 VECBOS 
prediction to fluctuate to the observed data is 1.8 x 10 -9 
(6cr for a Gaussian distribution). When the VECBOS 
calculation with Q 2 = is used, this probability is re- 
duced to 3.6 x 10 -6 (4.6cr). In the conservative case where 
Q 2 = M^r is used and the Et of each jet is increased by 
10% in the VECBOS program, the probability that the 
background is consistent with the W+ > 4-jet data is 
1.6 x 10" 4 {3.8a). 

To study the high-7^ events in the W+ > 4-jet sample, 
we perform binned maximum likelihood fits [11] of the 
data to a linear combination of the Ti. distributions pre- 
dicted by the VECBOS and the it HERWIG Monte Carlo 
calculations. Figure 3 shows the negative log-likelihoods 
for the two-component fits as a function of the top quark 
mass (Mtop)- After fitting the data points in this plot 
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FIG. 3. The least-squares fit of a cubic polynomial to the 
negative log-likelihood values from the two-component fits 
(see text), versus the top quark mass. The VECBOS Q 2 
scale is set to {Pt} 2 ■ The error bars reflect the statistical un- 
certainties of the fit due to finite Monte Carlo event samples. 
Also shown is the CDF mass result and error of Ref . [2]. 

to a cubic polynomial, we find the top quark mass to be 
180 ± 12(stat.)±^(syst.) GeV/c 2 . The systematic error 
reflects our uncertainties in the jet energy scale, the Q 2 
scale in VECBOS, and the level of initial state radiation 
predicted by the HERWIG it calculation. This top quark 
mass value is in excellent agreement with that found by 
the mass analysis of the 6-tagged W+ > 4-jet events [2]. 
The two estimates are correlated since the 19 events in 
Ref. [2] are a subset of the 99 events used here. Because of 
this, we only quote the Mtop estimate from this analy- 
sis to demonstrate its consistency with our earlier result 
of 176 ± 8 ± 10 GeV/c 2 , which remains our best mea- 
surement for Mtop- For a top mass of 180 GeV/c 2 , 
and setting Q 2 = (Pt) 2 in the VECBOS calculation, the 
fit yields 56 ± 10(stat.) ± 5(sys.) it events in the high- 
threshold signal sample of 99 events. The corresponding 



number with Q 



is 45 ± ll(stat.) ± 5(sys.). This 



can be compared with 34 ± 10(stat.) ± 5(sys.) events as 
extrapolated from the number of 6-tagged W+ > 3-jet 
events reported in [2]. These three estimates of the it 
production rate are consistent within errors. Table II 
summarizes the results of our analysis. 

We also perform a binned maximum likelihood fit to 
the low-threshold signal sample. The extracted top mass 
and the number of it events agree well with the results 
of the high-threshold fit, within statistical uncertainties. 

Figure 4 shows a single fit to the high-threshold signal 
sample assuming a top quark mass of 180 GeV/c 2 . The 
Q 2 scale was set to (P T ) 2 in the VECBOS prediction. 
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TABLE II. The fit results from the high-threshold signal 
sample for two choices of the Q 2 scale in VECBOS. The first 
error is statistical and the second error is systematic. The 
K.S. column is a Kolmogorov-Smirnov test confidence level. 



Q 2 


Mtop (GeV/c 2 ) 




K.S. (%) 


M 2 W 


180 ± 12+^ 


56 ± 10 ± 5 


56.5 


184 ± 15li| 


45 ± 11 ± 5 


63.2 



Good agreement between the best fit and the data is 
observed, as indicated by the Kolmogorov-Smirnov con- 
fidence levels shown in Table II. Also shown are the 23 
events tagged by the Silicon Vertex Detector (SVX) and 
the Soft Lepton Tagging (SLT) algorithms, out of which 
19 events had a good reconstructed fit to the it hypoth- 
esis in Ref. [2]. Most of the 6-tags lie in the region T~L 
> 270 GeV, which as the two-component fit shows, is 
dominated by it events. A study of background and sim- 
ulated it events shows that the dependence of the tagging 
efficiency on T~L is negligible. The clustering of 6-tagged 
events at large T~L provides additional evidence that the 
excess of events above the W+ jet background curve re- 
sults from it production. 

In summary, we compare the total transverse energy 
(Ti) distribution of W+ > 4-jet data with that expected 
from all known backgrounds and establish that they do 
not result from the same physical processes. We inter- 
pret the excess of high-7^ events observed in the data 
as the result of it production. The best fit to the 
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FIG. 4. The binned maximum likelihood fit of the 
high-threshold signal sample (solid line) to a linear combi- 
nation of the VECBOS W+ > 4-jet and HERWIG it predic- 
tions (dashed line), for Mtop = 180 GeV/c 2 . The dotted 
histogram is the it component of the fit. The Ti distribution 
for the SVX and SLT tagged events is also shown (shaded). 



data is a linear combination of W+ jet and it Monte 
Carlo calculations for a top mass of 180 GeV/c 2 . We also 
find a large fraction of events with at least one 6-tag in 
the W-enhanced high-7^ region, consistent with the stan- 
dard model prediction that one of the top quark decay 
products is the b quark. We demonstrate that a purely 
kinematic variable can be used to measure the top mass 
and the it production rate, and that Ti can be used in 
future analyses to discriminate against backgrounds to 
the top signal. 
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